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We develop an analytical approach based on a unitary transformation to investigate S — 1/2 antiferromag- 
netic Heisenberg chains coupled to phonons, and find a new quantum phase transition at zero temperature. 
Although the usual phase transition occurs depending on the strengths of the spin-spin and spin-phonon 
interactions, the phase transition in the present work is induced by a change in the geometrical structure of 
the spin-phonon interaction. The magnetic properties of the phase transition can be understood in relation 
to a phase transition that has already been found in the frustrated J1-J2 model. 

KEYWORDS: spin-phonon interaction, frustration, effective Hamiltonian, level spectroscopy, exact diagonaliza- 
tion, quantum phase transition, spin-Peierls compound 



1. Introduction 

Quantum spin systems with spin-phonon interactions 
have actively been studied as model systems for the re- 
search on spin-Peierls compounds, particularly the first 
inorganic spin-Peierls compound CuGcC>3. In the spin- 
Peierls compounds, a dimerized state with a spin gap 
due to a topological structure is realized at low tem- 
peratures. ^ It results from strong quantum fluctua- 
tions for quasi-one-dimensional systems. 2 ~ 4 ' In the one- 
dimensional systems, it has been theoretically elucidated 
that a spin-liquid state or a dimerized state is realized as 
the ground state depending on the strengths of the spin- 
spin and spin-phonon interactions, and the phase tran- 
sition between the states belongs to the same universal- 
ity class as that of the Berezinskii-Kosterlitz-Thouless- 
type 5 ' transition in the J1-J2 model. 6 ) 

The quantum spin systems with the spin-phonon in- 
teractions have theoretically been studied using effective 
spin Hamiltonians. 7-9 ' In the previous works, the effec- 
tive spin Hamiltonians are described by the J\-Ji model 
or the frustrated models with more long-range interac- 
tions. The Hamiltonian of the J\-Ji model is 

M = Ji Si ■ + J2 J]] Si ■ Si+2 , (1) 

i i 

where Si is the spin operator on site i. The values of J\ 
and J2 are positive and, therefore, the first and second 
terms on the right-hand side of eq. (1) are competitive in- 
teractions. In such a case, stable spin states are not easily 
found owing to the strong frustration. In a special case 
of J 2 = Ji/2, the ground state is exactly the product 
of spin-singlet states (Majumdar-Ghosh point). 10 ' The 
phase transition between the spin-liquid phase and the 
dimer phase in the J1-J2 model has been found by the 
level spectroscopy technique that considers the numerical 
data obtained by the exact diagonalization method. 11 ' In 
the quantum spin systems with the spin-phonon inter- 
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actions, besides these results, the following results have 
been obtained: both effects of phonons and frustration 
cooperatively work to form the gapped phase, 9 ' 12 ' and 
the separation of time scales of quantum fluctuations of 
spin and lattice degrees of freedom (the quantum nar- 
rowing effect) has been found by quantum Monte Cairo 
simulation. 13 ' 

In the present work, we have developed an analytical 
approach by which the phonon degrees of freedom are 
trace-outed and an effective Hamiltonian with frustration 
is yielded. Applying this method to a one-dimensional 
system with a more general spin-phonon interaction than 
usual cases of the previous works, we have found the 
ground-state phase transition between the spin-liquid 
phase and the dimer phase. In the usual cases, the spin- 
phonon interaction with difference coupling, 6 ~ 9, 14 ' 

?C = gJ2( b l + b ^ s * ■ s <+i - s < ■ s-i) > (2) 

i 

and/or that with local coupling, 9, 12, 13, 15, 

U% c = gJ2Q>l + h)Si-S i+1 , (3) 

i 

have been investigated, where b\ and bi are the phonon 
creation and annihilation operators at site i, respectively, 
and g is the coupling constant. The difference coupling 
was well studied before the discovery of CuGeC>3 and 
the local coupling was well studied in the research on 
CuGeOs. The magnetism of CuGe03 is due to Cu atoms, 
and the crystal structure shows that the spin-phonon in- 
teraction is described by 17 ' 

- E[-9 Cu ( 5 i+^- & i+i-^+i)+ff°' Go ( & i+^)]^-^+i . 

i 

(4) 

where Si denotes the spin of a Cu atom. As shown in 
Fig. 1, the spin-phonon interactions on Cu spins are in- 
duced by two types of displacement for Cu itself and O 
or Ge atoms. The term g Cu of eq. (4) indicates the spin- 
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Fig. 1. Illustration of lattice distortions for CuGe03. The spin- 
phonon interactions on Cu spins are induced by two types of dis- 
placement for Cu itself and O or Ge atoms. 



phonon interaction caused by the displacement of Cu 
atoms, and the term g°< Ge is caused by the displacement 
of or Ge atoms. Generally, although the displacements 
are written in terms of the different boson operators, we 
consider a simple mode: the Cu, O, and Ge atoms link 
and move, and the displacements are written in terms 
of one boson operator. In the research on CuGe03, the 
interaction without the term g Gu , i.e., the local coupling, 
is used as the relevant model because g°' Gc ^> g Gu . 17 ' 1S ) 

In the present work, we have studied the S — 1/2 spin- 
phonon interaction described by 



H sp = <?Y>] + h)(Si ■ S i+1 - ■ S^) 



(5) 



This interaction corresponds to interaction (4) through 
the relations 



gCu _j_ gO,Ge 



(6) 



and g = g Gn + g°> Gc . When .g°< Gc < g Gu and g°< Gc > 
<? Cu , (3 — 1 and 0, respectively. Therefore, this interaction 
includes both limits of the difference coupling (2) and 
the local coupling (3). In the present work, we have fo- 
cused on eq. (5) as the spin-phonon interaction and found 
the ground-state phase transition between the spin-liquid 
phase and the dimer phase by changing the value of j3. 
Although the phase transition is similar to that well- 
known in the system with usual interactions (2) and (3), 
it is a new phase transition in the sence that the driving 
force is /3 not g. 

The present paper is organized as follows. In §2, the 
analytical approach based on a unitary transformation 
is introduced and the effective Hamiltonians are ob- 
tained using the Ising and Heisenberg models coupled 
to phonons. In §3, numerical results for the S — 1/2 
antiferromagnetic Heisenberg chains coupled to phonons 
are presented. Section 4 is devoted to summary and dis- 
cussion. 

2. Effective Hamiltonian 

We consider the Hamiltonian described by 

h = h p + th B + gh sp , (7) 



i=l 
N 



X>1 + 6 *) V * > 



(8) 



(9) 



where N is the number of spin sites, g is the coupling con- 
stant of the spin-phonon interaction, and t is the positive 
exchange integral of the antiferromagnetic spin-spin in- 
teraction. The phonon frequency to is considered a unity, 
i.e., we assume an antiadiabatic case. The terms h s and vi 
have only to satisfy the conditions [h p , h s ] = [hp, Vi] = 0. 
We consider the following unitary transformation 9 ' for 
Hamiltonian (7): 



k=0 



p- ([<S , h p + th s ] k + k[S , /ispk-i) (10) 



fc! 

fe=0 

where S is an anti-Hermite operator of the order of g° 
and [A, B] k = [A, [A, ... [A, B] ]]. Wei/3e et al. 9 ~> de- 

s ^ 

k 

termined the form of S so that 'the amplitude of the 
first-order term in g of h operated to eigen states of the 
zeroth-order term in g is minimal.' In the present work, 
the specific form of S is determined so that 'the first- 
order term in g of h becomes zero': 



[S , hp + th s ] + h sp = 



(11) 



As a result, we could find the formal solution of eq. (11), 
S=[hp- th s , K] , (12) 

with 



/C — t 2n [h s , h S p]-2r, 



n=0 

Then, we obtain the transformed Hamiltonian 
' ' " J/. 



h = J2^{l-k)[S , h p + th s ] h 



k=0 



(13) 



(14) 



2.1 Ising model 

First, for a simple example, we consider the one- 
dimensional Ising model: 



h s = Y,(S?S? +1 + aS!S;, 



i+2J i 



and 



(15) 



(16) 



In this case, the transformed Hamiltonian is exactly ob- 
tained as 

h = h p +tJ2(S!S! +1 +a cS S^S^ 2 )-^g 2 (l+f3 2 ) , (17) 

i 

2 a 

where a c s — a + • The second term on the right-hand 
side of eq. (17) is the J%- J% Ising model. Its ground state 
is the Neel state for < a c g < 1/2 and the uudd state for 
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a c g > 1/2. 19 ) Thus, eq. (17) shows that the spin-phonon 
interaction acts cooperatively to form the uudd state for 
(3 > 0. When the ground state of original Hamiltonian (7) 
is expressed by |$q), we can obtain the expectation value 
of the lattice distortion ($ |&! + &i|$o) = 2g(l - fi)S 2 
for the Neel state. This result shows that the system is 
translated uniformly for (3 ^ 1. For the uudd state, on 
the other hand, ($ \b\ + k\$o) = (-l)*'2g(l + fi)S 2 , and 
the system is alternately distorted. 

2.2 Heisenberg model 

Next, we consider the one-dimensional S = 1/2 
Heisenberg model: 



Vi — SiSi+i — pSiSi-i 



(18) 



For simplicity, we set t = 0. The transformed Hamil- 
tonian up to the second order in g conserves the total 
phonon number. In this case, a vacuum state of phonon 
|0) p is an exact eigen state. Thus, the effective spin 
Hamiltonian h c g = p (0| h |0) p becomes 

,2 



leg 



9 -(l + f) 



(Si ■ s 



p 



i+1 



l + /3 : 



:Si-S 



16 



Ng 2 (l + P 2 ) • 



i+2) 



(19) 



We call it the 0(g 2 ) Hamiltonian. In the first term on 
the right-hand side of eq. (19), the ratio of the strength 
of the nearest-neighbor term to that of the next-nearest- 
neighbor term, which corresponds to J2/J1, does not de- 
pend on g. This is the J\-J% model with an effective 
frustration: 

P 



«cff 



1 + /3 2 



(20) 



For the local coupling with (3 — 0, a cS is zero and the 
0(g 2 ) Hamiltonian becomes the Heisenberg model with 
only the nearest-neighbor interaction. For the difference 
coupling with /3 — 1, a c s is 1/2 and the 0(g 2 ) Hamilto- 
nian describes the Majumdar- Ghosh point. 

The effective Hamiltonian up to the fourth order in g 
(we call it the 0(g A ) Hamiltonian) is 



toff = 



= Ji 



J\ J]] Si ■ Si+i + J2 2_. Si 

i i 
J3 ^ Si ■ Si + 3 



'i+2 



ri-Fl23 4 
T4(-F214 



F2-F213 
" F412 + 



f T3F312 
-F324 — -F423) 



where 



Fg- L s 2 s 3 



(Si ■ Si+s t )(Si-\ 



Jo = 



9i - j9i( l - a ff)(l + a ff) 



Ji 



9l ~ l^ti 1 ~ a cff)(l + «eff) 



(21) 
(22) 

N ,(23) 
(24) 



m 

1 

01. 

u? 




m 

1 

w 
'i- 




P 



Fig. 2. Comparisons between the ground-state energies of the 
original Hamiltonian and the effective Hamiltonians up to (a) sec- 
ond and (b) fourth orders in g for N = 10. 



IJ2 1 4 



<XeS 



r 2 
r 3 

r 4 



-(1 - a e fi)(l + a off ) 



:3i Q; off(2 + a e ff) , 



' 4 2 
rffi«off 



-^3i Q; cff(2 - 3a cff ) 
2S , i a cff(l - a c ff) , 



1 4 2 
-^5ia e ff 



(25) 
(26) 
(27) 
(28) 
(29) 
(30) 



and gi = g\/\ + (3 2 . The effective Hamiltonian could be 
described by gi and a e s only. Although the phonon vac- 
uum state |0) p is an exact eigen state of the 0(g 2 ) Hamil- 
tonian, this 0{g A ) Hamiltonian is an expectation value on 
|0) p . The third-nearest spin interaction and the multiple- 
spin exchange interactions arise from the fourth-order 
term in g. For /? = 0, however, the long-range inter- 
actions vanish and the 0(g 4 ) Hamiltonian becomes the 
J1-J2 model with the effective frustration: 

3 5 2 



n.- 



6g 2 



(31) 
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Fig. 3. Dependence of the ground-state energy per site for the 
original Hamiltonian on the cutoff parameter M for /3 = 1 and 
N = 10. The lines serve as visual guides. 



In the J\-Ji model, the phase transition between the 
spin-liquid phase and the dimer phase occurs at a c = 
0.2411. 20 ) Thus, for the 0( 5 4 ) Hamiltonian with (3 = 0, 

the phase transition occurs at g c = 



— / 8q c 



3+6q c 



0.66. 



3. Numerical Results 



To determine whether this effective-Hamiltonian 
method is useful, we show the j3 dependences of the 
differences between the ground-state energy per site 
Eorig/N for original Hamiltonian (7) and the E e s/N 
for effective Hamiltonians (19) and (21) in Fig. 2. 
The ground-state energies are calculated by the exact- 
diagonalization method with a periodic boundary condi- 
tion. In the original model, the total Hilbert space is writ- 
ten as the tensorial product space of spins and phonons. 
The dimension of the subspace associated to the phonons 
is infinite even for a finite chain system. Thus, we apply 
a truncation procedure retaining only basis states with 
the highest number of M phonons, where M is called 
the cutoff parameter. In the present work, we investi- 
gated the convergence of the ground-state energy for the 
truncation procedure. The dependence of E OI - lg /N on the 
cutoff parameter is shown in Fig. 3. The values of £" rig 
in Fig. 2 are obtained at M — 6. In Fig. 2(a), we show 
results for the 0(g 2 ) effective Hamiltonian. For g = 0.1, 
the energies -Eorig and E e s well agree in the entire range 
of < f3 < 2.0. When g increases, the difference in en- 
ergy increases for large f3 values. For g = 0.5, the energies 
disagree in the entire range of (3. In Fig. 2(b), we show 
results for the 0{g A ) effective Hamiltonian. In compari- 
son with the case of 0(g 2 ), the energies well agree in a 
wider range of /3. 

Next, we carry out level-spectroscopy analysis by using 
the numerical data obtained by the exact diagonaliza- 
tion of effective Hamiltonians (19) and (21). Using this 
method, Nomura and Okamoto 11 ) have analyzed low- 
excited states of the S = 1/2 antiferromagnetic XXZ 
spin chain with the next-nearest neighbor interactions 
and obtained the Berezinskii-Kosterlitz-Thouless-type 
phase-transition points from the energy-level crosses. As 
mentioned above, in the J1-J2 Heisenberg model with 
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P 

Fig. 4. p dependences of the lowest energy and the first-excited- 
state energy per site in the £;S? = and k = n subspace for 
g = 0.2 and N = 16. The square and circle represent the energies 
for the effective Hamiltonian up to the second and fourth orders in 
g, respectively. The lines serve as visual guides. 




0.005 0.01 
1//V 2 



0.015 



Fig. 5. Size dependence of /3 C (N) for various g's for the 0(</ 4 ) 
Hamiltonian. The solid lines are obtained by fitting the data of N = 
12, 16, and 20. The values at 1/7V 2 = are the phase-transition 
points in the thermodynamic limit. 



Hamiltonian (1), the phase-transition point is J2/J1 = 
0.2411. 20) Since our effective 0(g 2 ) Hamiltonian (19) is 
the same as the J1-J2 model Hamiltonian, the phase 
transition in the present work could also be evaluated by 
level-spectroscopy analysis. Although this is not trivial in 
the 0(<? 4 ) Hamiltonian with the multiple-spin-exchangc 
interactions, we could find an energy-level cross similar 
to that of the 0(g 2 ) Hamiltonian. Both level crosses are 
obtained using the lowest energy and the first-excited- 
state energy in the J2i = and k = it subspace, 
which correspond to the dimer and Neel excitations, re- 
spectively. 11 ) In Fig. 4, we show the (3 dependences of 
the lowest energy and the first-excited-state energy in 
the subspace. 

The phase-transition points in the thermodynamic 
limit are estimated by the least-squares fitting to the 
fitting function 20 ) 



A 
N 2 



(32) 
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Fig. 6. g dependence of the phase-transition point /3 C . The solid 
line corresponds to the transition point between the spin-liquid 
phase and the dimer phase for the 0(g 2 ) Hamiltonian. The circles 
represent the results for the 0(g 4 ) Hamiltonian. The dashed line 
serves as a visual guide. 



where A is a constant number. In Fig. 5, we show the 
system size dependence of (3 C for various <?'s for the 0(ff 4 ) 
Hamiltonian. The data of N — 12, 16, and 20 are fitted 
much better by the fitting function. 

In Fig. 6, we show the g dependence of the phase- 
transition point /3 C between the spin-liquid phase and 
the dimer phase obtained by level-spectroscopy analy- 
sis. For the 0(g 2 ) Hamiltonian, naturally, j3 c is con- 
stant. The value of 0.2571(1) well agrees with the solution 
/3 ~ 0.2570 of the equation a cB = 0.2411 = /3/(l + /3 2 ) 
obtained in the J\-Ji model. For large g values, we 
should consider the higher-order terms in g. For the 
0(g 4 ) Hamiltonian, f3 c decreases as g increases. This re- 
sult shows that the spin-phonon interaction acts cooper- 
atively to form the dimerized state. 

Thus far, we set t = for simplicity. If the contribution 
of the first-order term in t is taken into account under 
the condition t ~ g <C 1, the Hamiltonian th s = tyiASj- 
Si+i +aSi ■ Si+2) is added to the effective Hamiltonian. 
Namely, Ji+t and J2+ta replace J\ and J2, respectively. 
Then, the effective frustration becomes 



J 2 + ta J 2 



J2, 

+ {a-— ) 



t 



(33) 



Ji+t Jx x ~" Ji Ji+t 

The sign of the second term on the right-hand side of 
eq. (33) is positive for a > J2/J1 and negative for a < 
J2/ J\. Thus, the spin-spin interaction acts cooperatively 
to form the dimer phase for a > J 2 / J\ and the spin-liqud 
phase for a < J2I J\- 

We have also obtained results for the case in which we 
neglect the four-spin interactions in eq. (21). Since the 
results well agree with those for the non-neglecting case, 
the effect of the four-spin interactions is confirmed to be 
small. 

4. Summary and Discussion 

In the present work, we developed an analytical ap- 
proach based on a unitary transformation. Applying our 
method to the S = 1/2 antiferromagnetic Heisenberg 
chains coupled to phonons, we found a new quantum 



phase transition at zero temperature. Although the usual 
phase transition observed in the difference and local cou- 
pling cases occurs depending on the strengths of the spin- 
spin and spin-phonon interactions, the phase transition 
in the present work is induced by a change in the ge- 
ometrical structure of the spin-phonon interactions. We 
expect that it might be experimentally observed in new 
spin-Peierls compounds or as an effect of pressure. 

In the present work, we consider a nonadiabatic model 
and set t — for simplicity. Neutron scattering ex- 
periments have shown that neighboring spins within 
CuGeOs interact via high-cj phonons and the interac- 
tion is described by the nonadiabatic model with u ~ J, 
i.e., t ~ l. 15 ' 21 ) Furthermore, the strength of the spin- 
phonon interaction for CuGeC>3 has been estimated to 
be g ~ O.lt. 17 ^ Thus, the assumption of t — is not rea- 
sonable for CuGe03. Although the spin-spin interaction 
is an important factor for the research on CuGe03, 22 ) 
the result in the present work shows that the phase tran- 
sition occurs also in the system without the spin-spin 
interaction. The existence of the spin-spin interaction is 
not a necessary condition for the occurrence of the phase 
transition. 

The magnetic properties in the spin system coupled 
to phonons could be investigated from knowledge of the 
frustrated systems. Looking from a reverse viewpoint, 
we would be able to study the frustrated spin systems 
from knowledge of the spin systems coupled to phonons 
by quantum Monte Carlo simulation, the density matrix 
renormalization group method, and so on. 
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